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We explore an instantaneous decoherence correction (IDC) approach for the decoherence and energy relaxation 
in the quantum-classical dynamics of charge transport in organic semiconducting crystals. These effects, 
originating from environmental fluctuations, are essential ingredients of the carrier dynamics. The IDC is 
carried out by measurement-like operations in the adiabatic representation. While decoherence is inherent in 
the IDC, energy relaxation is taken into account by considering the detailed balance through the introduction 
of energy-dependent reweighing factors, which could be either Boltzmann (IDC-BM) or Miller-Abrahams 
(IDC-MA) type. For a non-diagonal electron-phonon coupling model, it is shown that the IDC tends to 
enhance diffusion while energy relaxation weakens this enhancement. As expected, both the IDC-BM and 
IDC-MA achieve a near-equilibrium distribution at finite temperatures in the diffusion process, while the 
Ehrenfest dynamics renders system tending to infinite temperature limit. The resulting energy relaxation 
times with the two kinds of factors lie in different regimes and exhibit different dependence on temperature, 
decoherence time and electron-phonon coupling strength, due to different dominant relaxation process. 


I. INTRODUCTION 

The advances of organic field-effect devices in re¬ 
cent years have inspired renewed research interest in 
the charge transport problems in organic semiconducting 
crystals.Recent studies focus on the intrinsic regime 
of charge transport, featured by observation of both 
band-like and hopping-like properties in the same de¬ 
vice under different operating conditions.Besides, the 
underlying electronic states are proved to be localized 
both theoretically'^ and experimentally.®dO R is realized 
that the coupling between carriers, which is denoted as 
electrons for simplicity throughout the paper, and lat¬ 
tice vibrations (phonons) plays an essential role in the 
charge transport process. Further, in many cases, the en¬ 
ergy scales of the electron-phonon coupling, temperature 
and intermolecular transfer integral overlap, which ren¬ 
ders the problem in the intermediated parameter regime, 
where the traditional pictures may not be effective.For 
example, Troisi and Orlandi proposed that the carriers 
in such materials are no longer localized by the self¬ 
trapping mechanism in polaron problems.Rather, it 
is due to dynamic disorder brought by phonons at fi¬ 
nite temperature. Many theoretical works have been 
contributed to this topic.Among them, the se¬ 
ries of studies using quantum-classical dynamics, which 
treat the electronic system quantum-mechanically and 
the phonons classically, are promising ones.^^ds They of¬ 
fer a time-dependent view of the transport process, en¬ 
abling the calculation of quantities such as the transient 
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conductivity.^'* Computationally, they are more efficient 
than the full quantum-mechanical methods and are con¬ 
venient to be tuned for realistic materials. For example, 
the detailed crystal and electronic structures extracted 
by experiments and first-principle calculations can be in¬ 
put as model parameters for multi-scale simulations.^** 
Besides, existing studies on chemical dynamics provide 
valuable resources for the improvement of those meth¬ 
ods. For example, Wang et al. successfully generalized 
the surface hopping approach to such systems,** reveal¬ 
ing an interesting crossover between band-like and hop¬ 
ping behavior. 

Despite the successes of the quantum-classical dynam¬ 
ics, there are still some inconsistencies in the resulting 
carrier dynamics. For example, in the pioneering work of 
Troisi et al. by Ehrenfest dynamics,*** the form of evo¬ 
lution of the electronic system is quantum-mechanically 
coherent. The electron wave functions keep expand¬ 
ing without a well-defined localization length, which 
does not coincide with experimental findings.**’***. This 
could be ascribed to the lack of explicit treatment of 
decoherence. ^®’'*‘*’'** In the charge transport process, the 
electronic system (electron) can be seen as the central 
system. In realistic situations, it would interact with 
many other species (environment), especially the domi¬ 
nant phonon degrees of freedom that are treated classi¬ 
cally here. Together, they form a system-environment 
architecture. One direct consequence for the system 
from coupling with environment is the destruction of 
phase coherence within the system, i.e., decoherence, 
which is essential for the consistency of quantum-classical 
dynamics.'*'* The incorporation of decoherence is indi¬ 
cated or proposed in several further studies. Fratini 
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et al. introduced a relaxation time approximation, in 
which the information of initial dynamics are used and 
the long-time dynamics are corrected.Yao et al. intro¬ 
duced the instantaneous decoherence correction (IDC) by 
measurement-like operations in the lattice site represen¬ 
tation, through which band-like temperature dependence 
of mobility and localized electronic states coexist.^' Also, 
in the surface hopping studies of Wang et al., decoher¬ 
ence corrections are included by collapsing the electronic 
states to the active states at the surface hopping events. 

Besides decoherence, another prominent process due to 
the system-environment coupling is the energy relaxation 
(dissipation) by which the system absorbs/dissipates 
energy from/to the environment and relaxes to the ther¬ 
mal equilibrium with the environmental temperature.^^ 
This process is essential for the interpretation of tran¬ 
sient experiments, such as the pump-probe ones.’^®’’^® 
However, there are signs that it is not treated properly 
in some forms of quantum-classical dynamics. For ex¬ 
ample, if the effect of phonon motion is taken to be a 
stochastic potential for the electronic system, all the adi¬ 
abatic states tend to be populated with equal probabil¬ 
ity in the long time limit.This means the effective 
temperature of the electronic system is infinitely large. 
Thus these methods are rigourous only when the ther¬ 
mal energy is much larger than the bandwidth. In the 
Ehrenfest dynamics the electronic system could influence 
the phonons by the mean-field Hellmann-Feynman force. 
However, it does not prevent the electronic system to 
deviate from the near-equilibrium distribution in the dif¬ 
fusion dynamics,which is also shown in the follow¬ 
ing. Besides, in the study using Kubo formalism with 
adiabatic approximation, the dc conductivity turns out 
to be zero.®^ The surface hopping approach is promis¬ 
ing in treating the energy relaxation through the frus¬ 
trated hops in the switching procedure.However, 
these hops also lead to inconsistency problems, the solu¬ 
tion of which is under active studies.®® Thus the problem 
of properly accounting for both decoherence and energy 
relaxation in the quantum-classical dynamics of the dy¬ 
namic disordered systems remains an open question. 

In this paper, we extend the IDC approach^' to the 
adiabatic representation to investigate the decoherence 
and energy relaxation processes. Similar approaches have 
been employed to the study of spin dynamics of excited 
states^^ and carrier dynamics of the Anderson model.**® 
The evolution of the system in short time regime (com¬ 
parable to the phonon frequency) are governed by the 
quantum-classical dynamics. The long time dynamics, 
which is recognized to be problematic,^® are modified 
by decoherence corrections. The IDC are carried out 
by measurement-like operations with different schemes. 
For energy relaxation, the detailed balance is considered 
by introducing the IDC with energy-dependent reweigh¬ 
ing factors, such as the Boltzmann (IDC-BM) and the 
Miller-Abrahams (IDC-MA) factors. The IDC with only 
destruction of phase coherence (IDC-DP) is also stud¬ 
ied for comparison. Based on the off-diagonal electron- 


phonon coupling model, the physical consequences of the 
IDC approach are explored. This paper is organized as 
follows. In Section H, the model and the IDC approach 
are presented in detail. In Section HI, the results of the 
IDC are shown, including both diffusion and electronic 
energy. The paper is summarized briefly in Section IV. 


II. MODEL AND METHOD 

In this section, we first present the off-diagonal 
electron-phonon coupling model and the equations of mo¬ 
tion. We then move on to introduce the IDC approach 
with a couple of schemes for energy relaxation, including 
that with Boltzmann and Miller-Abrahams factors. 


A. Model 


The Hamiltonian we consider in this work is composed 
of the electronic and the phonon part, which is H = 
Hei + ^Iph-*®’**** The electronic part is 

iJel = + Oi{Uj+i - Uj)]{c]cj+i -f c]_^iCj), (1) 

3 


where J is the transfer integral, a the electron-phonon 
coupling constant, Uj the displacement of phonon on the 
j-th site, cj (cj) the creation (annihilation) operators of 
electron. The Hamiltonian for the phonons is 


Ffph = 


E 

3 


\mu^j + 


( 2 ) 


where m is the effective mass of the phonon and ujq is 
the frequency. The electron state is described by the 
wave function |' 0 (t)) that is governed by the Schrodinger 
equation, which is 


(3) 

The electronic system back-reacts on the phonons 
through the mean-field Hellmann-Feynman force. The 
Langevin heat bath is included to account for the fluctu¬ 
ations of the phonons due to the thermal environment. 
The equations of motion for the phonons is 


2 

muj = —mLOMj --— 

dUj 


'ymuj + , 


( 4 ) 


where E^i = (*0(^)1 i?ei 1*0(0) is the expectation value of 
energy at time t, 7 the friction constant of the Langevin 
heat bath and the stochastic forces, satisfy the cor¬ 
relation function (0j(O'Ci(O)) = \/ 2 ^'rnkBTS{t)Sij, where 
ks is the Boltzmann constant and T is the tempera¬ 
ture. The initial displacements Mj( 0 ) and velocities Uj{ 0 ) 
are drawn from the Maxwell distributions with variance 
ksTjmujQ and ksTjm respectively. The initial state of 
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the electron is chosen among the adiabatic states (energy 
eigenstates) at t = 0 according to the Boltzmann distri¬ 
bution Pj, = exp{-E^/kBT)/J2^ey^p{-Eij,/kBT). The 
evolution is carried out by the 4th order Runge-Kutta 
method and the stochastic forces are incorporated by the 
method of Wang et al.^^’’*® The final result is averaged 
over enough realizations for the convergence of relevant 
physical quantities. 


B. Instantaneous Decoherence Correction 

We now move on to describe the instantaneous de¬ 
coherence corrections (IDC), which is carried out by 
repeated measurement-like operations on the electronic 
wave function in the adiabatic representation.^^ A de¬ 
coherence time td is introduced as the time interval be¬ 
tween successive IDC. It reflects the decoherence rate due 
to the system-environment interaction. This approach is 
shown to be effective to incorporate decoherence in the 
dynamics.Although discrete, it is shown in the site 
representation that the results resemble those by actually 
adding an attenuation term in the equation of motion of 
the non-diagonal elements of the density matrix.^' For 
the following results, a fixed td is used for convenience. 
Actually, more elaborate forms for td can be taken to 
reflect the situation in realistic materials. For example, 
proper distributions, such as the Poisson distribution,^^ 
can be used, together with specific models to determine 
td from other parameters. 

In detail, the dynamics are sliced into segments of time 
period td- Within each td, the carrier and phonons are 
evolved by the quantum-classical dynamics. After each 
td, a measurement-like operation is carried out for the 
electronic system by collapsing the electron state to a 
certain adiabatic state according to a chosen distribu¬ 
tion, which can be set by different schemes. Suppose the 
adiabatic states are denoted as |i?^) with energy and 
the distribution is P^, where /r is the index. A random 
number x is chosen uniformly in [0,1). The wave function 
is to be reset to that of \E,j) if 

V—1 V 

Y.p,<x<Y.p^- (5) 

/j.— 0 0 

For the distribution, the direct choice is that with the 
wave function in the adiabatic representation, which is 

p^ = I {E^ \m) l^ (6) 

where |' 0 (t)) is the wave function prior to the decoherence 
correction at time t. This scheme removes the phase 
relation among adiabatic states and is termed as the IDC 
with destruction of phase coherence (IDC-DP). However, 
as is shown below, this scheme only includes decoherence 
and energy relaxation is still not properly treated. 


C. Energy Relaxation 

For the purpose of energy relaxation in the IDC, we 
consider the detailed balance through the introduction 
of energy-dependent reweighing factors. We propose a 
couple of schemes in this work. The first one is to imple¬ 
ment a Boltzmann factor, which is 

-P™ = P^^ • exp(-P^/feBT)/CBM, (7) 

where Cbm is the normalization factor. This scheme 
could reflect the ability of the environment to recover 
the thermal distribution for the electronic system and 
is denoted as the IDC-BM in the following. It is moti¬ 
vated by Einstein’s theory of spontaneous emission of 
excited states,^® which is the result of quantum fluc¬ 
tuations of radiation fields. Further, motivated by the 
Miller-Abrahams formula,®*^ an alternative scheme is im¬ 
plemented. If the measurement is to collapse the state to 
an adiabatic state with higher energy, the probability is 
reduced by a factor from actually absorbing the energy 
from the environment. In this sense, the distribution is 

pMA _ / P^exp[-(P^ - E)/kBT]/CMA, Ef, > E, 

^ 1 PJCma, E^ < E, 

( 8 ) 

where Cma is a normalization constant. This scheme 
is termed as the IDC-MA in the following. Similar ex¬ 
pressions are widely used for the hopping rates in the 
simulation of charge transport in amorphous organic 
semiconductors.^ They are shown to capture the tran¬ 
sient mobility behavior'^^ and charge extraction tran¬ 
sients in organic solar cells,which depend on the correct 
treatment of energy relaxation. 

It is noted that with the IDC, the electronic states 
remain fairly localized. It is not necessary to solve the 
Hamiltonian of the whole lattice and a specific region 
is chosen instead, within which the electronic state is 
populated, which is more efficient. The region is deter¬ 
mined as follows. The index of the left/right bound¬ 
ary of the region is denoted as ji/r- The quantity 
Pl/r =J2f,Pt^\{jl/r is calculated and is ensured to 

be smaller than a critical value. Otherwise the region 
should be expanded. Actually, pi/^ reflect the expecta¬ 
tion value of populations on the boundary sites after the 
decoherence correction. In the following calculations, the 
critical value is chosen to be I 0 “®, which is small enough 
to not influence the final result. 


III. RESULTS 

For the following results, the parameters are taken to 
be typical for pentacene.^® The transfer integral J is 
0.037 eV, the electron-phonon coupling constant a = 3.3 
A“^, the phonon frequency Wq = 7.6 ps“^, the phonon 
effective mass m = 250 amu; the lattice constant a = 4 
A; the temperature T = 150 K. The friction constant 7 
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is taken to be 1 A lattice with 600 sites is taken 

to avoid the boundary effects. The decoherence time is 
taken to be td = IQh/ J w 180 fs. This set of parameters 
is taken in the following unless stated otherwise. 


A. DifFusion 

We first consider the diffusion with the IDC, which is 
reflected by the time-dependent averaged population of 
electron among lattice sites 

( 9 ) 

where th® coefficient of the electron wave function 

for the j-th site in the s-th realization. Ns is the to¬ 
tal number of realizations, which is taken to be at least 
10000. The mean squared displacement (MSD) is calcu¬ 
lated from the distribution Pj{t), which is 


MSD(t) = 




2 




( 10 ) 


The diffusion constant is the time derivative of MSD in 
the long time limit, which is 

I> = ll.„dMSDW 

2 t—>-oo dt 

The evolution time is taken to be at least 10 ps for the 
convergence of D. 

The typical MSD with the above IDC schemes are 
shown in Fig. 1(a), together with that from the Ehren- 
fest dynamics of Troisi et al.^''* The prominent effect of 
the IDC in adiabatic representation is that the diffusion 
is enhanced, in contrast to the suppression of diffusion 
with the IDC in site representation,^' which is a result of 
the quantum Zeno effect.The enhancement is reduced 
by inclusion of energy relaxation with the IDC-BM and 
IDC-MA schemes. For comparison, the results after av¬ 
eraging over 10 and 100 realizations with the IDC-BM 
are also shown. It can be seen that with 100 realizations, 
the result is already close to the final one. 

The enhancement of diffusion originates from the way 
the electron diffuses in the quantum-classical dynamics. 
Take the IDC-DP as an example. The average population 
on the j-th site before the decoherence correction is 

P^=\Y.{j\E;){E^\i^)\\ ( 12 ) 

The decoherence correction destroys the phase coherence 
between different adiabatic states and the corresponding 
population after the decoherence correction is 

P^ = Y,\{j\E^)\^\{E^\^)\\ (13) 
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FIG. 1. (a) Evolution of the mean squared displacement 

(MSD) with different instantaneous decoherence correction 
(IDC) schemes applied to an off-diagonal electron-phonon 
coupling model, in comparison with the result of Troisi’s 
method with Ehrenfest dynamics (green dot-dot-dash): IDC 
with Boltzmann factor (IDC-BM, red solid), IDC with Miller- 
Abrahams factor (IDC-MA, blue dash) and IDC with destruc¬ 
tion of phase coherence (IDC-DP, purple dot-dash). All re¬ 
sults are got by averaging over more than 10000 realizations. 
The results with the IDC-BM with 10 and 100 realizations 
are also shown (grey solid), (b) Temperature dependence 
of diffusion constant with the above IDC schemes: IDC-BM 
(red square), IDC-MA (blue circle), IDC-DP (purple lower- 
triangle) and Ehrenfest dynamics (green upper triangle). 


Comparing with P^, the missing terms in P^ are the in¬ 
terference ones between adiabatic states. In the Ehren¬ 
fest dynamics, the wave function starts to change by the 
mixing of adiabatic states due to phonon motion. If 
the wave function is expanded in the evolved adiabatic 
states, the interference terms tend to be constructive for 
those sites that are initially populated and destructive 
for those sites that are initially unpopulated. With fur¬ 
ther evolution, this kind of phase relation is disrupted by 
the stochastic phonon motion and population emerges 
on those initially unpopulated sites. The IDC helps this 
process and thus enhances diffusion. A specific exam¬ 
ple of this enhancement is shown in the supplementary 
material.’'* 

We further calculate the temperature dependence of 
diffusion constants. The results are shown in Fig. 1(b). 
It can be seen that all the IDC schemes, like the Ehrenfest 
dynamics, give a band-like dependence of diffusion con¬ 
stant. In the lower temperature regime around 100 K, the 
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FIG. 2. (a) Evolution of the electronic energy with the 
Ehrenfest dynamics starting from the adiabatic states with 
lowest (red dash) and highest (blue dot-dash) energies, and 
energies that are closest to ±1.5J, ±J, ±0.5J (grey solid). 
The result starting from random adiabatic states with equal 
probability (green solid) is also shown, (b) Evolution of Egi 
for the IDC-BM (red solid), IDC-MA (blue dash) and IDC- 
DP (purple dot-dash) schemes starting from adiabatic states 
with the lowest and the highest energies. The average value 
of the Boltzmann distribution is shown by the arrow on the 
left. 


diffusion constant with the IDC are substantially larger 
than the ones given by the Ehrenfest dynamics. With 
increasing temperature, the diffusion constants from the 
IDC-DP scheme remain higher than the Ehrenfest ones, 
while those from the IDC-BM and IDC-MA schemes 
gradually decrease and nearly coincide with the Ehrenfest 
ones. Besides, the diffusion constants with the IDC show 
sharper dependence on temperature. Beyond 150 K, they 
begin to deviate from a power-law behavior. This can be 
related to the corresponding deviation of the localization 
length of adiabatic states. The temperature dependence 
of the averaged localization length is provided in the sup¬ 
plementary material.^'* 


B. Energy Distribution of Electronic States 

We now turn our attention to the electronic energy 
Ee\ = (V^(t)l He\ \'ip{t)) for the energy relaxation process. 
Firstly, the evolution of Eei in the Ehrenfest dynamics 
are shown in Fig. 2(a). The initial condition of the elec¬ 
tronic system is chosen to be adiabatic states with both 
the highest and the lowest energies, and those closest to 
±1.5J, ±J, ±0.5J. They all tend to relax to energy 


values near the center of the density of states {E = 0). 
To check the observation, the evolution of Eg\ starting 
from random adiabatic states with equal probability is 
also shown. An initial small decrease away from zero is 
observed, which persists beyond 10 ps. The small devia¬ 
tion is caused by the polaron effect due to the Hellmann- 
Feynman force, which is a feedback from the electronic 
system to the phonons. It is confirmed by the evolution 
of averaged displacements and phonon potential energies 
that are shown in the supplementary material.^’'* These 
results show that the near-equilibrium Boltzmann distri¬ 
bution is violated in the Ehrenfest dynamics. All the 
adiabatic states tend to be populated with equal prob¬ 
ability, which corresponds to a large effective tempera¬ 
ture. This violation is from the inability of the coherent 
quantum dynamics to distinguish between the processes 
that go upward and downward in energy. Similar prob¬ 
lems are also pointed out in the studies of the spin-boson 
dynamics by stochastic Schodinger equation with real¬ 
valued stochastic terms^^ and the Anderson model with 
Haken-Strobl-Reineker method.^® 

We then study the evolution of the electronic energy 
with the IDC schemes starting from adiabatic states with 
both the highest and the lowest energies, which are shown 
in Fig. 2(b). For the IDC-DP scheme, the final energy 
tends to a finite negative value different from the value 
of the Boltzmann distribution, which is indicated by the 
arrow on the left. This negative value has the same origin 
with the negative dip in the Ehrenfest dynamics as dis¬ 
cussed above. Here, the effect is made more pronounced 
by the decoherence corrections to adiabatic states. Fur¬ 
ther, it is clear that for the IDC-BM scheme, the elec¬ 
tronic energies starting from both the lowest and highest 
adiabatic states tend to the same limit, which is only 
slightly lower (~ 0.06 J) than the Boltzmann value. For 
the IDC-MA scheme, the energies also relax to a final 
value, which is slightly higher than the Boltzmann one. 
Besides, the IDC-BM scheme gives a much faster relax¬ 
ation process than the IDC-MA. In all, both schemes help 
recover the near-equilibrium Boltzmann distribution in 
the dynamics, which is not observed in the Ehrenfest dy¬ 
namics and the IDC-DP scheme. This property is present 
in the temperature regime considered in this paper, which 
is shown in the supplementary material. 

We further calculated the population distributions 
P{E) of the the electronic states with respect to energy 
after evolving the system for 10 ps, which are shown in 
Fig. 3(a). In the calculation of P{E), the contribution 
from each realization is broadened by a Lorentzian line- 
shape with width P, which is 

P{E) = CpJ2 \{K\*‘)f r; . (14) 

S,/i '' 

where Cp is the normalization constant to ensure 
f P{E)dE = 1. Here P is chosen to be 0.05J. The den¬ 
sity of states D{E) can be defined similarly. A quantity 
f{E) is further defined to depict the probability distri¬ 
bution of the electronic system as f{E) = P{E)/D{E), 
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FIG. 3. (a) Population distribution P{E) of the electronic 
states with respect to energy E after evolution of 10 ps start¬ 
ing from adiabatic states with the Boltzmann distribution: 
IDC-BM (red solid), IDC-MA (blue dash), IDC-DP (purple 
dot-dash) and Ehrenfest dynamics (green dot-dot-dash). The 
result of the Boltzmann distribution (black triple-dot-dash) 
and the density of states D{E) (black dot) are also shown 
for comparison, (b) Corresponding probability distribution 
defined as f{E) = P{E)/D{E). 


which is shown in Fig. 3(b). In both figures, the energy 
is from — 2.5J to — J, in which the population is pro¬ 
nounced. A plot over the full energy regime is provided 
in the supplementary material for clarity. For compar¬ 
ison, the results of the Boltzmann distribution are also 
shown. The deviation of the Boltzmann result from lin¬ 
ear behavior around E = —2.2J in Fig. 3(b) is an ar¬ 
tifact from the fast variation of D{E). In Fig. 3(a), 
D{E) is also shown as the case when all the adiabatic 
states are populated with equal probability. It can be 
seen that the results of the Ehrenfest dynamics tend to 
D{E) and the results of the IDC-DP scheme also de¬ 
viate from the near-equilibrium Boltzmann distribution. 
In contrast, the results of both the IDC-BM and IDC- 
MA are not far away from equilibrium ones. Besides, 
in the lower energy regime below the peak of the pop¬ 
ulation, the IDC-BM gives a distribution that is closer 
to the Boltzmann one; in the higher energy regime, the 
correspondence is better with the IDC-MA scheme. 


C. Energy Relaxation Time 

Furthermore, we note that an energy relaxation time 
E can be extracted from the the evolution of electronic 
energy, which is a quantity commonly used for describ- 



r(K) 



td (fs) 



FIG. 4. (a) Dependence of the inverse energy relaxation time 
on temperature T for the IDG-BM (red square) and the 
IDC-MA (blue circle) schemes. The decoherence time is cho¬ 
sen to be td ~ 180 fs. (b) Dependence of with decoherence 
time td- The temperature is fixed at T = 150 K. (c) Depen¬ 
dence of t~^ with electron-phonon couplings a, with oq = 3.3 

A-i. 


ing the associated transient properties.^"* For this pur¬ 
pose, we calculate the evolution of Eei from randomly 
chosen adiabatic states with equal probability. The evo¬ 
lution can be well fitted by an exponential damping 
function exp(—t/t^) in the temperature regions con¬ 
sidered. The dependence of the inverse energy relax¬ 
ation time on temperature, decoherence time and 
electron-phonon couplings are shown in Fig. 4, where 
large means faster energy relaxation. It can be seen 
that the behaviors of with the IDC-BM and IDC- 
BM schemes are very different. For the temperature de¬ 
pendence, from the IDC-BM scheme decreases with 
temperature while that from the IDC-MA scheme in¬ 
creases. For the decoherence-time dependence, from 
the IDC-BM scheme decreases with td while that from 
the IDC-MA scheme is insensitive. Both schemes give in- 
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verse that increase with increasing electron-phonon 
coupling a. These differences come from the different 
dominant energy relaxation processes. In the IDC-BM, 
energy relaxation is dominated by the decoherence cor¬ 
rection. With lower temperature, the Boltzmann factor 
~ exp{—E/kBT) has greater ability to redistribute the 
population toward lower adiabatic states. For smaller 
td, the redistribution is more frequent. In both cases, 
the relaxation is faster. In contrast, for the IDC-MA 
scheme, energy relaxation is through the energy fluctu¬ 
ation of adiabatic states. The decoherence corrections 
bring about the imbalance between the processes that 
go upward and downward in energy. With lower tem¬ 
perature, the fluctuation is slower, leading to slower en¬ 
ergy relaxation. Besides, different td do not influence 
the fluctuation substantially, thus leading to the insen¬ 
sitivity. For both schemes, increasing electron-phonon 
coupling leads to faster mixing of adiabatic states and 
faster energy relaxation. Besides, it is noted that tr with 
the IDC-MA (I ps) are generally ne order of magnitude 
larger than the implemented decoherence time td (100 fs), 
which follow the same relationship revealed in the stud¬ 
ies of some elementary models, such as the spin-boson 
model.®® In contrast, tr with the IDC-BM are compara¬ 
ble to td^ which may correspond to a faster relaxation 
mechanism.®® 


IV. SUMMARY 

In summary, we have studied the decoherence and en¬ 
ergy relaxation, which are essential processes induced 
by the (quantum) system-environment interaction in the 
quantum-classical dynamics of charge transport in or¬ 
ganic semiconducting crystals. They can be incorpo¬ 
rated by the IDC with either the Boltzmann or the 
Miller-Abrahams factors. We find that, while decoher¬ 
ence enhances diffusion, the energy relaxation weakens 
this enhancement. With both the IDC-BM and IDC-MA, 
the distributions of electronic states tend to the near¬ 
equilibrium one, which is a sign of proper treatment of 
energy relaxation. The energy relaxation time from the 
two schemes show different behaviors in its dependence 
on temperature, decoherence time and electron-phonon 
coupling strength, which result from different relaxation 
processes. Furthermore, it can be shown that the deco¬ 
herence and energy relaxation are crucial for observing 
a direct current response with quantum-classical dynam¬ 
ics, which make a direct calculation of mobility under a 
finite external electric field possible.®® 
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